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Abstract 

We study non-topological solitons, so called Q-balls, which carry a non- vanishing Noether 
charge and arise as lump solutions of self-interacting complex scalar field models. Explicit 
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examples of new axially symmetric non-spinning Q-ball solutions that have not been studied 
so far are constructed numerically. These solutions can be interpreted as angular excitations 
of the fundamental Q-balls and are related to the spherical harmonics. Correspondingly, 
they have higher energy and their energy densities possess two local maxima on the positive 
z-axis. 

We also study two Q-balls interacting via a potential term in 3 + 1 dimensions and 
construct examples of stationary solitonic-like objects in (3 + l)-dimensional flat space- 
time that consist of two interacting global scalar fields. We concentrate on configurations 
composed of one spinning and one non-spinning Q-ball and study the parameter-dependence 
of the energy and charges of the configuration. 
, In addition, we present numerical evidence that for fixed values of the coupling constants 

two different types of 2-Q-ball solutions exist: solutions with defined parity, but also solutions 
^3 i which are asymmetric with respect to reflexion through the x-y-plane. 

o\ 

1 Introduction 

i> 

Solitons play an important role in many areas of physics. As classical solutions of non-linear 
field theories, they are localised structures with finite energy, which are globally regular. In 
general, one can distinguish topological and non-topological solitons. While topological solitons 
[X] possess a conserved quantity, the topological charge, that stems (in most cases) from the 
spontaneous symmetry breaking of the theory, non-topological solitons [21 [3] have a conserved 
Noether charge that results from a symmetry of the Lagrangian. The standard example of 
non-topological solitons are Q-balls [3], which are solutions of theories with self- interacting 
complex scalar fields. These objects are stationary with an explicitely time-dependent phase. 
The conserved Noether charge Q is then related to the global phase invariance of the theory and 
is directly proportional to the frequency. Q can e.g. be interpreted as particle number [2]. 
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While in standard scalar field theories, it was shown that a non-renormalisable <]? -potential 
is necessary [5], super symmetric extensions of the Standard Model (SM) also possess Q-ball 
solutions [6]. In the latter case, several scalar fields interact via complicated potentials. It was 
shown that cubic interaction terms that result from Yukawa couplings in the superpotential and 
supersymmetry breaking terms lead to the existence of Q-balls with non-vanishing baryon or 
lepton number or electric charge. These supersymmetric Q-balls have been considered recently 
as possible candidates for baryonic dark matter |7j and their astrophysical implications have 
been discussed [8]. 

Two interacting scalar fields are also interesting from another point of view. Up until now, 
the number of explicit examples of stationary solitonic-like solutions that involve two interacting 
global scalar fields is small. An important example are superconducting strings, which are axially 
symmetric in 2 + 1 dimensions extended trivially into the z-direction |9j. Axially symmetric 
generalisations in 3 + 1 dimensions, so-called vortons, have been constructed in [10]. Note that 
all these solutions have been constructed in models which have a renormalisable $ 4 -potential. 

Here, we study two interacting scalar fields in 3+1 dimensions and construct explicit examples 
of stationary solitonic-like axially symmetric solutions consisting of two global scalar fields. 
While vortons possess one scalar field with an unbroken U(l) symmetry (the condensate field) 
and a scalar field whose U(l) is spontaneously broken (the string field), we here consider two 
scalar fields with unbroken U(l) symmetries. One can thus see our model as the limit of vanishing 
vacuum expectation value for the second scalar field. Then, stationary solitonic like objects can 
be constructed explicitely. Note that the model in [10] contains a renormalisable ^-potential, 
while we need a non-renormalisable $ 6 -potential here. However, as stated in [10] , the explicit 
construction of vortons was done using also a non-renormalisable potential which contains an 
interaction term of the form <&\<&\. 

Q-ball solutions in 3+1 dimensions have been first studied in detail in [5]. It was realised that 
next to non-spinning Q-balls, which are spherically symmetric, spinning solutions exist. These 
are axially symmetric with energy density of toroidal shape and angular momentum J = kQ, 
where Q is the Noether charge of the solution and k £ Z corresponds to the winding around the 
z-axis. Approximated solutions of the non-linear partial differential equations were constructed 
in [5] by means of a truncated series in the spherical harmonics to describe the angular part of 
the solutions. The full partial differential equation was solved numerically in [TT]. It was also 
realised in [3] that in each A;-sector, parity-even (P = +1) and parity-odd (P = —1) solutions 
exist. Parity-even and parity-odd refers to the fact that the solution is symmetric and anti- 
symmetric, respectively with respect to a reflection through the x-y-plane, i.e. under 9 —* it — 9. 

These two types of solutions are closely related to the fact that the angular part of the so- 
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lutions constructed in OUT] is connected to the spherical harmonic Yq(0,ip) for the spherically 
symmetric Q-ball, to the spherical harmonic Y-±(6,tp) for the spinning parity even (P = +1) 
solution and to the spherical harmonic Y^d^fp) for the parity odd (P = —1) solution, respec- 
tively. Radially excited solutions of the spherically symmetric, non-spinning solution were also 
obtained. These solutions are still spherically symmetric but the scalar field develops one or 
several nodes for r e]0, oo[. In relation to the apparent connection of the angular part of the 
known solutions to the spherical harmonics, it is natural to investigate whether "^-angular ex- 
citations" of the Q-balls exist in correspondence to the whole family of spherical harmonics 
Yj?(9,tp), —L < k < L. This can further be motivated by the fact that, in the small field limit 
where a linear approximation can be used, the field equation describing the Q-ball becomes a 
standard harmonic equation that can be solved by separation of variables and whose fundamen- 
tal solutions are given in terms of spherical harmonics for the angular part. Of course, it has to 
be checked whether this correspondence, expected from the linear limit, still holds for the full, 
i.e. non-linear equation. 

In the present paper, we present strong numerical arguments that new angularly excited 
solutions of the non-linear field equations exist and that the correspondence between angular 
excitations of the Q-balls and spherical harmonics indeed holds. In addition to the solutions 
corresponding to Y£ and Y^~ x for k = 1,2,3 presented in [5] we have constructed solutions 
with angular dependence and symmetries corresponding to the spherical harmonics Y® and Y" 2 °. 
These solutions are non-spinning but constitute axially symmetric excitations with respect to 
the angular coordinate 0. As expected, these new solutions have higher energies and charges 
than the spherically symmetric solutions and we would thus expect them to be unstable. These 
solutions thus complete the already known spectrum of Q-ball solutions and show that not only 
radial excitations of fundamental soliton solutions, but also angular excitations exist. 

We also study two interacting Q-balls and put the emphasis on the interaction between a 
non-spinning and a spinning Q-ball. In particular, we investigate the dependence of the energy 
and the charges of the solution on the interaction parameter and the frequencies, respectively. 

Next to parity-even and parity-odd solutions, we also construct solutions that have no defined 
parity with respect to reflexion through the x-y-plane. 

The explicit construction of solutions with two interacting complex scalar fields is surely of 
interest for the astrophysical implications of such objects, especially for the construction of such 
objects in supersymmetric theories. Moreover, it adds to the spectrum of soliton solutions that 
e.g. possess no definite parity. 

The differential equations describing both excited as well as interacting Q-balls are non-linear 
partial differential equations, which -to our knowledge- cannot be solved analytically. We thus 
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solve these equations numerically using an appropriate PDE solver |12j . 

Our paper is organised as follows: in Section 2, we discuss the model and give the equations 
and boundary conditions. In Section 3, we discuss the new Q-ball solutions for k = 0, while in 
Section 4, we present our results for two interacting Q-balls. Section 5 contains our conclusions. 



2 The model 

In the following, we study a scalar field model in 3 + 1 dimensions describing two Q-balls 
interacting via a potential term. The Lagrangian reads: 

c = -d^d^i + ia^a^a - v($i, ^) (i) 

where both <3?i and $2 are complex scalar fields. The potential reads: 

2 

[/($!, $ 2 ) = N^l 6 - Pi\$i\ 4 + l^i\ 2 ) + A|$i| 2 |$ 2 | 2 (2) 
i=i 

where aj, 7$, i = 1, 2 are the standard potential parameters for each Q-ball, while A denotes 
the interaction parameter. 

In [5] it was argued that a <I> 6 -potential is necessary in order to have classical Q-ball solutions. 
This is still necessary for the model we have defined here, since we want $1 = and $2 = to 
be a local minimum of the potential. A pure $ 4 -potential which is bounded from below wouldn't 
fulfill these criteria. 

The Lagrangian ([T]) is invariant under the two global U(l) transformations 

$1 -> $ie iai , $ 2 $ 2 e ia2 (3) 

which can be applied separately or together. As such the total conserved Noether current ff tat \ , 
fi = 0,1,2,3, associated to these symmetries is just the sum of the two individually conserved 
currents j± and j% with 

tftot) = ft + # = -« ind^i - ®id»$$) - i (*5S^* a - * a ^*5) . (4) 

with = 0, = and d^ tot) = 0. 

The total Noether charge Q(tot) °f the system is then the sum of the two individual Noether 
charges Q\ and Q2: 

Q(tot) = Qi + Q2 = - J fi<$x - J j° 2 d 3 x (5) 

Finally, the energy-momentum tensor reads: 

2 

T^ u = {d^id^* + d u ^d^l) - g^C (6) 
i=i 



2.1 Ansatz 



We choose as Ansatz for the fields in spherical coordinates: 

$i(t, r, 9, ip) = e^+^Ur, 0) , i = 1, 2 (7) 

where the and the k{ are constants. Since we require $i(<p) = &i(<p + 27r), i = 1,2, we 
have that ki £ Z. It was moreover demonstrated in j5j [TT] that Q-balls exist only in a specific 
parameter range uj m i n < cu < uj max and that the charge Q tends to infinity when either u — > w m in 
or u —* u> max . We discuss the limits in the 2 Q-ball system in the following section. 
The Noether charges of the solution then read: 

/fir roc 
\<S>i\ 2 (fx = A-KUJi I I r 2 sin 9 dr d9 $ , i = l,2 (8) 
Jo Jo 

while the energy is given by the volume integral of the tt-component of the energy-momentum 



tensor: 

TT OO 

E= / Tnn d 3 x = 2ir 







r 2 sin9 dr d9^l Uifa + (^) 2 + 



i=l 



+ -^L + "i|0i| 6 -^i| 4 + 7iW 2 ) +A|<Ai| 2 |02 12 
r z sin 9 



0) 



where the prime and dot denote the derivative with respect to r and 9, respectively. 

For ki 7^ 0, the solutions have non-vanishing angular momentum that is quantised. The total 
angular momentum J is the sum of the angular momenta of the two individual Q-balls: 

J = j T 0ip d 3 x = J 1 + J 2 = k x Qx + k 2 Q 2 (10) 

We will thus in the following refer to solutions with ki = as non-spinning and to solutions with 
ki as spinning. 

The Euler-Lagrange equations read: 

2 11 k 2 

<t>" + -4 + ^<f>i + ^ cot 94>i - ! 2 6i + uj 2 4>i = 3ai<t>! ~ + li4>i + Xfctl ( n ) 

r r z r z r z sin 9 

with % = 1, 2 and k ^ i. 

The boundary conditions, which result from requirements of regularity, finiteness of the 
energy and the symmetry of the solution, are: 

d r <f>i(r = 0,9) = , do<f>i(r = oo,0)=0, d g <f>i(r, 9 = 0, tt) = , i = l,2. (12) 

for non-spinning solutions with ki = and 

<k(r = 0, 9) = , <pi{r = oo, 6) = , &(r, = 0, tt) = , i = 1, 2 . (13) 

for spinning solutions fcj 7^ 0. 
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2.2 Bounds on U\ and ui 2 in the 2-Q-ball system 

In j5j [11] the bounds on the frequency oj have been discussed in the case of one Q-ball. Here, we 
note that these bounds have to be modified if one considers two interacting Q-balls. The set of 
equations (llip can be interpreted as the mechanical equations describing the frictional motion 
of a particle in two dimensions. The effective potential in this case reads: 

V{<h,<h) = \{<4<ft + - \u{<k,<h) (14) 

Q-ball solutions exist provided the configuration (0i = 0,02 = 0) corresponds to a local max- 
imum of the effective potential and provided the effective potential has positive values in any 
radial direction from the origin in the <j>i — ^2— plane. This leads to non-trivial bounds for the 
parameters uj\ and 0J2- 

The former condition leads to the requirement that 

w l < ^l,max =71 , ^2 < ^Imax = 72 • (15) 

The latter condition leads to a more complicated domain of existence in the u;i-L<J2-pla n e. To 
describe this condition, we introduce the polar decomposition of <p\ and 02 as follows: 

0i=pcosx , 02 = psinx • (16) 

where < x < 2-7T and < p < oo. 

The condition on the frequencies uj\ and 002 then read: 

u\ cos 2 x + ^l sin 2 x > (^1 cos 2 x + ^l sin 2 x)min = ram p [U(p,x)/p 2 } , Vx (17) 

In the particular case that we have studied throughout this paper, namely a\ = ct2 = 1, 
Pi = P2 = 2 and 71 = 72 = 1.1 this inequality takes the form: 

lo\ cos 2 x ~t~ ^2 sin 2 x > [— 5A 2 cos 4 x sin 4 x + 20A cos 2 xsin 2 x(cos x + sin x) 

.0 o 9 (\ 9 

+ 2(cos x + sin x + Hcos x sm X+Hsin xcos x 

— 20 cos 4 x sin 4 x)] / (cos 4 x + sin 4 \ — cos 2 x sin 2 x) (18) 

For x = titt/2, n = 0,1,2,..., we recover the results of the one Q-ball system discussed in 
O E]. For all other values of x, the limiting values for lo\ and u>2 will depend on the value of 
the interaction coupling A. E.g. for 0i = 02, i-e. x = ^/A, we find : 

uj\+u? 2 > 1/5 + A- 1/8A 2 (19) 

Thus, for small A, the lower bound on the value of tof + a; 2 will be larger than in the non- 
interacting limit. 
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3 New non-spinning Q-ball solutions for 012 = 02 = 72 = A = 



In order to be able to understand the structure of a system of two Q-balls, we have reconsidered 
the one Q-ball system. We set all quantities with index "2" to zero in the following and omit 
the index "1" for the remaining quantities. 

In this section, we would like to point out that more than the previously in the literature 
discussed solutions exist. 

For this, we first consider the equation for one Q-ball with vanishing potential. This reads: 

2 11 k 2 

<f>" + ~(f>' + -^<t>+ -j cot 9<j)- „ , 2 6 + uj 2 4> = (20) 
r r z r z r z sin 

Although the solutions of the above equation are well known, it will be useful for the following 
to recall their properties. Using the standard separation of variables, the solutions read: 

*(r,M« JL+l/ f r) Yt(eM (2D 

where J denotes the Bessel function, while Y£ are the standard spherical harmonics with 
-L < k < L. 

One may hope that solutions of the full equations with the discrete symmetries corresponding 
to the ones of the spherical harmonics will exist. Of course, the non-linear potential interaction 
will deform the radial part of the solutions of the linear equation in a highly non-trivial manner. 

The solutions of the full equation constructed so far for k = have been spherically sym- 
metric. With the above arguments, axially symmetric solutions should equally exist with an 
angular dependence of the form Y®, e.g. for L = 1, the angular dependence should be of the 
form cos 9. In the following, we will denote the solutions of the full non-linear equations with 
angular symmetries corresponding to the symmetries of the spherical harmonic Yj? by <j)\. 

3.1 Numerical results 

The partial differential equation has been solved numerically subject to the boundary conditions 
fl 121) or (|13p using the finite difference solver FIDISOL [12] . We have mapped the infinite interval 
of the r coordinate [0 : 00] to the finite compact interval [0 : 1] using the new coordinate 
z := r/(r + 1). We have typically used grid sizes of 150 points in r-direction and 50 points in 
8 direction. The solutions have relative errors of 10~ 3 or smaller. Throughout this section, we 
choose a\ = a = 1, 0\ = (3 = 2, 71 = 7 = 1.1. 

In FigH] (left), we show the profile of a new solution that we obtained for k = and 
uj = 0.8. This solution looks like a deformation of the spherical harmonic Y® with the appropriate 
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Figure 1: The profile of the function (f>® is shown for u = 0.8, a = 1, j3 = 2, 7 = 1.1 (left). The 
corresponding energy density Too is also given (right). 

symmetry with respect to 9 = tt/2 and is clearly axially symmetric. In particular it fulfills 
0?(r,7r/2)=O. 

The field (j)\{r,9) is maximal at a finite distance from the origin on the positive z-axis. 
Moreover, the configuration is anti-symmetric under reflexion through the x — y— plane, i.e. 
under 9 —* tt — 9. Thus the solution is parity-odd: P = —1. Note that we have only plotted the 
function for 9 6 [0 : tt/2], but that we have verified the symmetry of the solution. 

We also present the corresponding energy density Too m Fig{T] (right). It shows that the 
density of the solution is mainly concentrated within two small "balls" situated around the 
positive z-axis (at z ~ 2.4 and z ~ 7.6) and separated by a minimum (at z ~ 5). The position 
of this minimum coincides with the maximum of the scalar field (4>\)max ~ <Ai(5,0) ~ 1.2. It 
can be checked that this value corresponds roughly to a local minimum of the potential while 
the partial derivatives are evidently small in this region, explaining the occurrence of a minimal 
value of the energy density at (x,y,z) (0,0,5). Of course, due to the anti-symmetry of the 
solution this pattern is equally given on the negative z-axis. 

The classical energy and charge of this new solution is higher than that of the spherically 
symmetric k = solution (see Table 1 and Table 2 below) , however lower than that of the k = 1 
spinning Q-ball. 

In order to investigate further our idea of constructing new solutions as deformations of the 
spherical harmonics, we have also investigated solutions with higher value of L and we managed 

8 











\\ . 












\ 










\ 




\\ 


2 


*,° \\ 


3 

1 : cos( 8) 


\\ 


■ V 


2 : (3 cos( e) 2 - 1)/4 


3 : (5cos( 0) 3 -3cos( 0) )/2 

i i i i 


I I 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 





Fi gure 2: (f>® (solid), 2 (short dashed) and (f>® (long dashed) are shown as functions of 6 for 
fixed value of r. We have chosen r ~ 5 for (f)\, r ~ 2 for 2 and r ~ 6 for 0g. The corresponding 
spherical harmonics oc cos(0), Y 2 ° ^ 3cos(#) 2 — 1 and Y® oc 5cos(6>) 3 — 3cos# (with an 
appropriate normalisation) are also shown. 

to construct solutions 2 and (j)® corresponding in their angular symmetries to those of the 
spherical harmonics Y 2 ° °^ 3 cos 2 6 — 1 and Y® oc 5 cos 3 6 — 3 cos 6, respectively. 

In Fig|2j we plot (f>®, (fy and 4>z as functions of 6 for a fixed value of r together with the 
corresponding spherical harmonics Y®, Y® and Y®. Here, we have chosen r ~ 5 for ifi®, r ~ 2 
for 2 and r ~ 6 for <f>^. The first thing to notice is that the symmetries of the solutions , 
4>2 an d 03 with respect to reflection at 6 = ir/2 are exactly equal to those of the corresponding 
spherical harmonics. The actual solutions are, of course, deformed with respect to the spherical 
harmonics, but the correspondence is apparent. E.g. the solution 2 has <900 2 (r, 7r /2) = (in 
contrast to the solution 0° which has (^^(r, tt/2) = 0). We don't show the energy density of 2 
and 03 here, since it resembles that shown in Fig. 1. 

We believe that the correspondence also holds for higher spherical harmonics. 

Since we have presented strong numerical evidence that the correspondence with the spherical 
harmonics holds, it is justified to label the different solutions of the field equation by means of 
the quantum numbers of the corresponding spherical harmonic, i.e. by L and k referring to 
with L, k integers and —L<k<L. Needless to say that the numerical construction becomes 
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Figure 3: The profile of the function (/>!] for u = 0.8, a = 1, f3 = 2, 7 = 1.1 



more involved when the difference L — \k\ increases. Adopting these notations and fixing the 
potential according to a = 1, /3 = 2, 7 = 1.1, we find for the solutions corresponding to w = 0.8 
and to = 0.84 the values for the energy E and charge per frequency — given in Table 1 and Table 
2, respectively. 

The first three solutions L = 0,1,2 in this list are static (i.e. non spinning) while the 
last, <p\, is stationary (i.e. spinning). For all the solutions we constructed, the energy of the 
non-spinning solutions is lower that the energy of the spinning ones. 
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symmetry 
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Table 1: The energy E, the charge per frequency Q/u>, the parity P and the symmetry of the 
first few Q-ball solutions is given for uj = 0.8. 
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Table 2: The energy E, the charge per frequency Q/u>, the parity P and the symmetry of the 
first few Q-ball solutions is given for u> = 0.84. 

4 Interacting Q-balls 

Since in supersymmetric extensions of the Standard Model, Q-balls exist that result from the 
interaction of several scalar fields, we here investigate the interaction of two classical Q-balls as 
toy-model for these systems. 

For two spherically symmetric Q-balls (k\ = k<i = 0) in interaction, the 2-Q-ball solution is 
still spherically symmetric and the domain of existence in the cji-cj2-plane can be determined by 
using the reasoning given in Section 2.2. Here, we put the emphasis on solutions where the two 
Q-balls have different symmetries and study the effect of the direct interaction parameterized 
by the coupling constant A. 

We believe that a particularly interesting case is the interaction between a spherically sym- 
metric, non-spinning Q-ball (k\ = 0) and a spinning Q-ball (&2 = 0). We have thus restricted 
our analysis to this case and set ki = and k-z = 1 in the following. 

Note that we will index all quantities related to the spherical Q-ball in the following with 
"1", while all quantities related to the axially symmetric Q-ball will be indexed with "2". 

For later use, we define the "binding energy" of the solution according to 

AE = E — Ek 1= o — Ek 2 =\ ■ (22) 

It represents the difference between the energy E of the 2-Q-ball configuration and the sum of the 
energies of the two single (i.e. non-interacting) Q-balls Ek 1= o, -Etc 2 =i with the same frequency. 
We expect that those solutions which have AE < to be stable, while those with AE > would 
be unstable. 

4.1 Numerical results 

We have solved the two coupled partial differential equations using the solver FIDISOL [12] for 
several values of u>i, u>2 and A and fixing a% = a.2 = 1, /?i = 02 = 2 and 71 = 72 = 1.1. As 
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Figure 4: The energy density Too of the 2-Q-ball solution consisting of a spherically symmetric, 
non-spinning Q-ball (k\ = 0) and a spinning Q-ball (&2 = 1) is shown for uj\ = L02 = 0.8, 
oc% = Q!2 = I, Pi = P2 = 2, 71 = 72 = 1-1 and for three different values of A = 0, 1, —0.5. 

starting profiles, we have used the corresponding non-interacting Q-ball solutions. For A = 0, 
these solve the two decoupled partial differential equations. We have then slowly increases the 
parameter A to obtain the interacting solutions. 
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4.1.1 UJ i = U>2 



In order to understand the influence of the interaction parameter A, we show the energy density 
Too f° r u>i = UJ2 = 0.8 and three different values of A in FigJH For A = 0, the two Q-balls are non- 
interacting and the energy density is just a simple superposition of the energy densities of the 
two individual Q-balls. For A 7- the Q-balls interact. For A > 0, it is energetically favourable 
to have the two Q-bahV cores in different regions of space. As seen in FigJH for A = 1, the 
spinning Q-ball seems to be "pushed away" from the non-spinning, spherically symmetric one. 
For A < 0, it is energetically favourable to have two Q-balls sitting "on top of each other" . This 
is shown in FigUJfor A = —0.5, where the two Q-balls seem to be localised at the same place. 

We have also studied the dependence of the energy E, the binding energy AE and the 
two charges Q\ and Q2 on the interaction parameter A. The results are shown in Figj5] for 
U\ = o»2 = 0.8. All quantities increase with the increase of A, specifically it is evident that the 
2-Q-ball configuration is energetically more stable for A < than for A > 0. Specifically, we 
would thus expect the solution to be stable for A < and unstable for A > 0. 

Following our discussion in Section 2.2, we have also studied the dependence of the energy E 
and of the charges Q\ and Q2 on the frequencies u% and u>2- Our results for uj\ = U2 are shown 
in Figj6]for A = —0.5, and 0.5, respectively. 
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Figure 5: The quantities E, Qi, Q2 and AE are shown as functions of the interaction parameter 
A with = L02 = 0.8, ai = 02 = 1, /?i = /?2 = 2, 71 = 72 = 1.1 

As expected, the energy E for a given frequency uj\ = LO2 is higher (resp. lower) than in the 
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Figure 6: The energy E (left) and the charges Q\ and Q2 (right) are shown as functions of the 
frequency uj\ = U2 for A = —0.5, 0, 0.5. 

non-interacting case for positive (resp. negative) values of A. 

As before, we find that the solutions exist in a given interval of the frequency: u;i jm j n (A) < 
^ < ^>i,maxW (and equally for u>2 since u>i = U2). We have determined the bounds on uj\ and 
L02 in Section 2.2 for two spherically symmetric Q-balls. Here, we would expect that these values 
change slightly since we have a system of one spherically symmetric and one axially symmetric 
Q-ball. However, we see that the qualitative results are similar here. We observe that for A > 0, 
the values of the energy E and of the charges Q±, Q2 diverge at u\ = u)\^ m i n and uo\ = u>i )max . 
Following the discussion of Section 2.2 we find that the maximal value of uj\ is independent of 
A. This can be clearly seen in Fig. where the energy E and the charges Q\ and Q2 diverge 
at u)\ — &i,max ~ 1.035 for all three values of A. Note that this maximal value is only slightly 
lower than the bound given in Section 2.2: uj\ max = 1.1. The reason why the bound is not equal 
is that here we are dealing with an axially symmetric solution interacting with a spherically 
symmetric one. Analytic arguments of the type done in Section 2.2 are, however, only possible 
if the Euler-Lagrange equations are ordinary differential equations, i.e. only in the case where 
the solutions are spherically symmetric. So, it is not surprising that the analytic values differ 
from the numerical ones. 

On the other hand, the minimal value of uj\ is A-dependent. This can be seen in Fig. We 
have given our results only for u > 0.6 in this figure since the construction of solutions becomes 
increasingly difficult for uj < 0.6. However, it can be clearly seen that the energy E and Q2 
diverge at different values of u)\ — ^i,mm* In agreement with Section 2.2., we find that u)\^fii^i is 
increasing for increasing (and small) A. 
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For A < the behaviour at the lower bound of uj\ changes. We observe that Qi correspond- 
ing to the spherically symmetric field <j)\ decreases when io\ decreases. The analysis of the profile 
of the solution reveals that the field <p\ deviates only slightly from the spherically symmetric 
configuration for frequencies close to u)\^ max - However, it gets more and more deformed in the 
equatorial plane when uj\ decreases. At the same time, the field </>2 increases in the equatorial 
plane. This phenomenon is illustrated in FigJ7]for A = —0.5, to = 0.6 and u = uJi tTnax , respec- 
tively. In this figure, the fields (pi and <p2 as well as the energy density Too are shown as function 
of r for two angles 9 = and 9 = tt/2, respectively. 

4.1.2 Solutions with lo\ ^ uj2 

We have also constructed 2-Q-ball solutions for u% ^ uj2- The energy density Too of a 2-Q-ball 
solution corresponding to ui = 0.65 and L02 = 1 is shown in Figj8] for four different values of 
A. The result is qualitatively similar to the case uji = u>2- This figure however suggests very 
clearly that for A < the &2 = 1 Q-ball has tendency to disappear from the 2-Q-ball system. 
For instance the maximal value of the 4>2 field, \(j)2,max\, decreases for decreasing A. 

We have also studied the dependence of the solution's conserved quantities on L02 = o;i/0.65 
for A = ±0.5. The dependence of the energy E and the charges Qi, Q2 is shown in Figj9j These 
results strongly suggest that for A < and in the region of the parameter space chosen, the field 
(p2 corresponding to the k2 = 1 Q-ball tends uniformly to zero for a critical value of u>2 = £<4 
such that Q2 — > for ui2 — ► t<4 . Only the field (\)\ remains non-trivial when u>2 < c<4 • This 
effect can also be observed in FigiEJ where the solution for A = —0.5 has nearly lost all its axially 
symmetric character. 

We observe the inverse phenomenon for ui = 002 with a constant c > 0. We don't present 
our detailed results here since they are qualitatively equivalent to the case discussed above. We 
find that Q\ — > for uj\ — > u^f r \ Thus, the spherically symmetric solution disappears from the 
system, while $2 remains non-trivial for u\ < Lu^f r \ 

Apparently, while in the case uj\ = UJ2 and A < 0, the charge Q\ associated to the spherical 
Q-ball tends to zero for u\ — ► uj^f r \ it is the charge Qi, i = 1,2 of the Q-ball with the higher 
frequency that tends to zero for u>, — > u;j cr \ i = 1,2 when oj\ ^ 102 and A < 0. Note that nothing 
similar is observed when A > 0. 

While 1-Q-ball solutions known so far are always either parity-even or parity-odd with respect 
to 9 — » 7r — 9, we have constructed several examples of 2-Q-ball solutions that do not have a 
defined parity. One such solution is shown in FigllOl (lower part) together with a parity-even 
solution (upper part). These solutions exist for exactly the same values of the coupling constants. 
Both functions <j>\, 4>2 are clearly neither parity-even nor parity-odd and the field 4>2 possesses in 
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Figure 7: The profiles of <pi, <p2 and the energy density Too are shown for 9 = and 9 = ir/2, 
respectively. Here A = —0.5. The upper figure is for to = 0.6, the lower for u> = uJ\, m ax ~ 1-035. 



addition nodes in the radial direction. This solution is thus an asymmetric, radially excited 2- 
Q-ball solution. As expected, we observe that this asymmetric solution has much higher energy 
and charges than the corresponding parity-even solution. 

The investigation of solutions of this type and their eventual bifurcation into branches of 
solutions with defined parity is currently underway. 
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Figure 8: The energy density Too of the 2-Q-ball solution consisting of a spherically symmetric, 
non-spinning Q-ball (k\ = 0) and a spinning Q-ball (fo = 1) is shown for uj\ = 0.65, 1V2 = 1, 
ot% = «2 = lj 0i = 02 = 2, 71 = 72 = 1-1 and for four different values of A = 0.5, 0, —0.2 and 
-0.5. 

5 Concluding remarks 

In this paper, we have presented numerical evidence that Q-balls solutions admit several types 
of excitations labelled by integers. So far, it was known that the static, spherically symmetric 
solution is the "ground state" of a series of radially excited solutions. Families of spinning 
solutions are also known, they are axially symmetric and can be labelled according to the 
winding k around the axis of symmetry. Here we present evidence that excitation with respect 
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Figure 9: The energy E and charges Q\, Q2 are shown as functions of 002 for lo\ = 0.65u;2 and 
A = ±0.5. 

to 9 can be constructed as well. Generally, the previous results and the present analysis suggest 
that families of elementary solutions of the field equations exist and are labelled by n, L, k, where 
n refers to the number of nodes in radial direction, while L, k refer to the "quantum numbers" 
related to the spherical harmonics. At the moment, the only analytic argument we have for this 
property is its analogy to the linearized version (i.e. small field limit) of the equation where this 
result holds true by standard harmonic analysis. It is likely that the qualitative properties of 
the solutions exist also in the case of the full non-linear equations. 

We have also studied a system of two interacting Q-balls and have constructed several exam- 
ples of axially symmetric, stationary solutions that carry conserved currents and charges. We 
observe that the 2-Q-ball solutions exist in a finite range of the frequency 0Ji, m in < U{ < 0Ji^ max , 
i = 1,2, where 0Ji^ max is independent of the interaction coupling, while uJi, m in does dependent of 
the interaction coupling in a highly non-trivial manner. We find that the charges Qi, i = 1, 2 of 
the 2 Q-balls in interaction tend to infinity when Ui — ► oJi^max or Ui — > Ui^min as long as A > 0. 
For A < 0, however, we observe that the charges Qi associated to the Q-ball with the higher 
frequency uoi tends to zero for Ui — > < uJi^ max . For u>i im m < &i < uf r ^ only the remaining 
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Figure 10: The contour plots for <pi and 4>2 of a parity-even 2-Q-ball solution (upper part) and of 
an asymmetric 2-Q-ball solution (lower part) are shown for A = —0.5, u% = 0.585 and uji = 0.9 

field (j)j, j ^ i is non-zero. 

In a future publication, we intend to construct solutions with the more realistic potential 
available from supersymmetry [6] and put the emphasis on the possibility of constructing Q-balls 
and their excited and/or spinning versions with potentials involving only quartic terms in the 
scalar fields. 
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